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Abstract 

Origin of ultra high energy cosmic rays is an unsolved problem. 
Several proposals such as Z-burst, decay of super massive matter, 
susy particles as a primary, neutrino as a primary in extra dimen- 
sion models exist in the literature which try to address this issue. 
Many of these proposals solve the problem of propagation of cosmic 
rays over cosmological distances by introducing new physics. However 
these do not explain the origin of such high energy cosmic rays. The 
possible astrophysics sites, such as active galactic nuclei, are highly 
constrained. Here we determine whether these cosmic rays originated 
from the decay of some exotic objects, such as primordial black holes 
(PBHs), present in the early universe. In contrast to the usual Top 
Down scenario we do not assume that this exotic object necessarily 
has to decay in our astrophysical neighbourhood since we assume a 
beyond the standard model scenario, where the propagation problem 
is absent. We consider the standard 4-dimension PBHs as well as the 
brane world PBHs. We find that in both cases it is unable to produce 
the observed ultra high energy cosmic ray flux. 
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1 Introduction 



The observation of cosmic rays with energies in excess of 10 20 eV present 
a major challenge to astro-particle physics. Due to the presence of cosmic 
microwave background radiation it is predicted that cosmic rays with ener- 
gies above 10 20 eV will not be observed due to the GZK cut-off £Q. The 
presence of GZK violating events implies new physics unless the source of 
these events lie within our astrophysical neighbourhood. Some interesting 
possibilities include topological defects [3 E] , primordial black holes [I] or 
super heavy particles decaying within a distance of about 100 Mpc in 
order to evade the GZK bound. However most of such topdown scenarios 
are severely constrained by existing data. There also exist many propos- 
als which solve the cosmological propagation problem by introducing new 
physics. Examples include violations of lorentz invariance El, existence 
of susy particles |7j, existence of magnetic monopoles jH], Z-burst [3110], a 
strongly interacting neutrino at ultra high energies(UHE) [TT] etc. Many of 
these proposals are also severely constrained by existing experiments and will 
be further tested by future planned experiments. The strongly interacting 
neutrino proposal, for example, will be ruled out by non-observation of UHE 
neutrinos in experiments by the year 2006 10 . 

The only possible conventional astrophysical sources for UHE neutrinos 
are Active Galactic Nuclei (AGN) and Gamma Ray Bursts (GRB). These can 
be considered as possible sources of ultra high energy cosmic rays (UHECR) 
only if we assume that the propagation problem is solved by some new physics 
proposal. It is clearly important to consider alternate sources of UHECR, 
even if they are located at cosmological distances. Here we consider Primor- 
dial black holes decaying in the early universe. These objects are interesting 
because, depending on their masses, they can survive till today. UHECRs 
from Primordial black holes(PBH) in our astrophysical neighbourhood have 
been studied in Ref. jl]. 

In this paper we consider the production of UHE protons and neutrinos 
from PBHs decaying today and also PBHs which decay in early epoch of the 
cosmological evolution of the universe. We calculate UHE fluxes in standard 
4D PBHs as well as in 5D braneworld PBHs. In the next two sections we 
review the 4D PBHs and 5D brane world PBHs. 



2 



2 Standard 4D Primordial Black Holes 



It is known that black holes would have formed in very early universe through 
density fluctuations [12] • These fluctuations may either be primordial or may 
be formed spontaneously at any epoch. PBHs might have formed at any 
spontaneous symmetry breaking epoch through the collision of bubbles |13 J 
or through the collapse of cosmic strings ^3]. If one assumes the production 
of black hole with mass of the order of horizon mass at some time t in the 
evolution of universe then its mass [15] , 

M BH {t) « tl „ 10^ (*) g. (1) 



G V10- 23 s, 

Such a black hole with mass of order 10 15 g would be evaporating now. Masses 
less than 10 15 g would have evaporated by now. For example, black holes with 
masses between 10 10 g and 10 13 g would have completed their evaporation 
between 10 3 and 10 12 sec. 

The emission rate of spin 1/2 particles from a black hole is given by [TE] . 



d 2 iV r 1/2 (£,T) 



dtdE exp (J.) + i 



(2) 



per particle degree of freedom. Here T is the temperature of the hole and 
Ti/ 2 (E,T) is the absorption coefficient. The above spectrum deviates from 
the black body spectrum due to the energy dependence of T 1 / 2 (E,T). In 
the limit Ej (kT) ^> 1, the spectrum approaches that of a black body with a 
temperature JTH] 

ig 



T w 1.06 x 10 13 , (3) 

[Ml v ; 

where M is the mass of the black hole. For relativistic particles Ti/ 2 {E,T) 

is given by [THJ HZ] , 

27 E 2 



A black hole loses mass at the rate fUJ E] 

dM a(M) 



(5) 



dt M 2 ' 

where a(M) depends on the degrees of freedom of the emitted particles and 
increases with temperature. For standard model 

a(M) w lO 26 ^ 3 ^- 1 (6) 
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above top quark production threshhold. From eq. |3] and eq. 03 we find jl] 



dU = 1.5 x 10- 15 — *- , (7) 

* 

where t* = t/(lsec) and T* = T/(l EeV). Particles with energies above 
1 EeV will be produced instantly when the temperature of the black hole 
is such that kT > 1 EeV. Characteristic time for the production of EeV 
energy particles is of the order of 10~ 18 sec Similarly characteristic time 
required to produce particles with planck energy 10 9 EeV is of the order of 
10~ 43 sec. Hence duration of the emission of these high energy particles by 
PBHs is unimportant compared to the evolution time of the universe. After 
integrating eq. the mass of the black hole, at any time t, is approximately 

given by psmansi, 

where Mi, to, M* are respectively the initial mass of the PBHs, age of the 
universe and mass of the PBHs evaporating today. The Standard Model 
provides a lower bound on the high energy value of a(M), above top pro- 
duction threshold. One can assume same value of a(M) to hold at higher 
energies. In next section we review the properties of brane world PBHs. 



3 Braneworld Primordial Black Holes 

Braneworld cosmological models provide an interesting alternative to the 
standard cosmology. In this scenario PBHs can be formed in very early uni- 
verse by density perturbation [T^]. In the well known RS2 model our universe 
is a positive tension brane in a 5d bulk with s 1 jZ-i type of compactification. 
It is shown in ref. |2H] that in this model PBHs can be easily produced in 
the radiation dominated universe through spherical collapse. 
The Friedmann equation in RS2 model is ^J] 

tt<2 8vr / p 2 \ A 4 k 

where H is the Hubble constant, p, p are respectively the energy density and 
pressure of the fluid, a is the scale factor on the brane, M 4 is the effective 
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4D planck mass, A 4 the 4D cosmological constant and pkk the dark radiation 
density. The energy conservation equation on the brane leads to 



p + 3H(p + p) = . (10) 

As usual k = —1,0, 1 for open, flat and closed brane universe. From nucle- 
osynthesis constraints one finds that the pkk term in equation is negligible. 
The energy density and scale factor then becomes (THj 

3M| 1 

P = ; r 11 

H 32tt t(t + t c ) y ' 

1/4 

(12) 



and 




respectively, where to is any nonzero time, t c is the transition time, t c — |, 
and I is AdS curvature radius. The time t c separates nonconventional cos- 
mology from the standard cosmology. For time t ^> t c we recover stan- 
dard cosmology. For time t <C t c , which we call high energy regime, the 
temperature-time relation modifies to [19] 

T /45X 1 / 4 (IY 1,A ft\-^ 



T 4 V87r 3 / " c \UJ \U 

where T 4 , t 4 , and Z 4 are 4-D Planck temperature, Planck time and Planck 
length respectively. At transition time t c , transition temperature is [T9| |2*T] 

T c = 3xl0 18 f^j GeV. (14) 

Since current experiments probe gravity to a length scale of 0.2 mm and 
constrain the size of extra dimension to be / < 10 31 / 4 , we find that the 
minimum value of T r is 10 3 GeV. 



3.1 Evaporation rate of Brane World PBHs 

In ref. authors have calculated a mass-lifetime relation for black holes 
formed on the brane due to collapse of matter on the brane. Effect of ac- 
cretion on the lifetime and mass of brane world PBHs are studied in ref. 
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[2H I2H1 12I| • We do not consider the effect of accretion in this paper. For a 
review on brane world PBH, see ref. (2H|. If the size of the black holes r <C / 
then their geometry is described by 5D schwarchild black holes. These black 
holes will emit hawking radiation into the brane as well as to the bulk. In 
this approximation, radius and temperature of the black hole are given by 

1X9 > 

i\ 1 ' 2 fM^i 2 



and 




= I H [ W J h (15) 



T bh = — (16) 
2nr 

respectively. Mass loss rate of these black holes is given by [TH] . 

dM 16tt, 

dt ~ 

where 



gT , (17) 



1 9C(5) 

9 — TT^gbrane + QO , gbulk ■ (18) 

loll oZtt 

In our case g& rane dominates and most of energy goes to the brane. If we 
consider only the Standard Model degrees of freedom, gbrane = 100. From eq. 
[T71we can derive the lifetime t evap of a black hole of initial mass M. Lifetime 

tevap is 

^"' ibs) ■ (19) 

eq. El gives a relation between time and temperature of the BH as 

.OO^' 1 dT* 
512ttA ~fl> 



.poor 1 dr. 

dt * = ^ , (2°) 



where i* = t/lsec, T* = T/l EeV and A = j^. The parameter A basically 
determines the length I at which 5D BHs dominate the dynamics and can 
atmost take value 10 31 [22] • This value comes from upper limit on the size 
of extra dimension constrained by sub-mm gravity experiments. From mass- 
lifetime relation we get a range of / over which PBH acts as a 5 dimensional 
black hole. The minimum value of AdS radius (THj is 

ft \ V 3 

Imin = ~ 9 1/3 (=) h (21) 
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and the maximum value of AdS radius l max = 10 31 Z 4 . This implies that the 
5-dimensional PBHs lie in the range ^J] 

// \ -V2 /f s 1/2 

M mm = f/ 2 r-^J (^J M 4 (22) 

to 

M max = ^/ 3 f%^V /3 M 4 . (23) 



v t 4 / 

Similarly temperature of the PBHs ranges from 

T mi „ = J^-g- 1/4 f ^1 7%^) ~ 1/4 T 4 (24) 



327T y Imin J V t 4 

to 

We have listed a range of mass, size and temperature of brane world PBHs 
corresponding to their evaporation in Table. 



4 Mass distribution of black hole 

The mass distribution functions for 4D PBHs formed by the collapse of den- 
sity perturbations have been derived in ref. ^3]. Here we assume that the 
mass function is dominated by a particular value, as given in j2Hj. In this 
case the distribution of the PBHs present throughout the evolution of the 
universe from their time of formation can be expressed as, 

n(M) = N(l + zf . (26) 

where z is the redshift at time t. Here we have neglected the evaporation 
of PBHs from their formation time to their time of evaporation. We assume 
similar mass function for 5D brane world PBHs. 



5 Flux Calculation for 4D PBHs 

We are interested in the flux of particles, specifically neutrinos and protons, 
at ultra high energies from PBHs. The neutrinos may be emitted directly 
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by the PBH, which we refer to as the direct flux. Alternatively they may be 
emitted by hadrons and other particles in the PBH spectrum. We refer to 
this as the indirect flux. 



5.1 Direct Neutrino Flux 

Let f(E v ,T) represent the direct neutrino flux at energy E u given in eq. 121 
then the diffuse flux per unit area today is: 

_ i xl . 5xl0 - r r d M^^f(K,T^, (27) 



dE M 4tt J Zmin JkT„(i+z) (kT*) 4 (1 + z) 2 dz 

where E v = E u0 (l + z), z is the redshift at the time of emission, z max cor- 
responds to maximum redshift from which particles of energies 100 EeV can 
reach us and E u0 is energy of neutrino at redshift z min . We take z max < 10 7 , 
which corresponds to the redshift at which particles of initial energy 10 19 GeV 
will be observed at energy of 100 EeV. 



5.2 Indirect Neutrino flux 

The indirect neutrino flux is obtained dominantly from the decay of hadrons. 
The hadrons are formed by the fragmentation of quarks and gluons. The 
hadronic flux is generally dominated by pions which form almost 97 % of the 
flux. The remaining 3 % is mostly nucleons. The fragmentation function for 
quarks into hadrons is given by 

where £ p = Y{\ + (c /Y)^ - c /Y) and 2a 2 = (|rj) 1/2 with Y = 
ln(^),b = m^nz, Co = _jg_ ari d a = if- + f^. The symbols N c 
and np refer to the number of colors and number of flavours respectively. 
For iV c = 3 and rip = 6, the parameter b is equal to 7. We fix the normal- 
ization constant by equating the multiplicity of corresponding hadrons to 
their experimental value at Q = \fs = 91 GeV and A e // = 200 MeV. For 
pions we find Ah = 4.89. 

In order to calculate indirect neutrino flux we consider the following pro- 
cesses. 
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Figure 1: Direct Neutrino flux today from 4d PBHs evaporating at redshift 
z = O(solid), z = 1000(short dashed) and z = 10 6 (long dashed). Sim- 
ilarly indirect neutrino flux today from 4d PBHs evaporating at redshift 
z = O(dottted), z = 1000 (small spaced dots) and z = 10 6 (large spaced 
dots). 



1. the decays of and pions. 

2. the fragmentation of quarks into pions and then through the following 
channel 7r — > ji — > v. 

3. the decays of evaporated W-bosons through the following channel W — > 
e + v and W — > ji — > z/. 

While calculating the indirect neutrino flux we have incorporated the 
fragmentation functions by modifying the expression f(E u ,T). The final ex- 
pression for f(E v ,T) in this case is given in the Appendix of this paper. 
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Figure 2: Proton flux today from 4d PBHs evaporating at redshift z = 
(solid), z = 1000 (dotted) and z = 10 6 (small spaced dots). Similarly 
Proton flux today from 5d BWPBHs evaporating at redshift z = (large 
spaced dots), z = 1000(long dashed) and z = 10 6 (short dashed). 



5.3 Proton Flux 

The expression for the proton flux is similarly given by 



dN p 



dE p0 An 



1x1.5x10-- f™ [ tT - rfJt 1 ^HE„T)^dz, 
t7r J Zmin JkT u (i+z) Je p {kT*y (1 + zy dz dE p 

(29) 

where E p = E p0 (l + z), z is the redshift at the time of emission, ^ = .03^-, 
x = -j*- and E p q is the energy of proton at redshift z min . 
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6 Flux Calculation for 5D brane world PBHs 



A 5D black hole emits particles with energy in the range (E, E + dE) at a 
rate 

= 1 E2 (30) 

dtdE 47r3 T 2 exp(f)+l 1 ' 

per particle degree of freedom. Here T is the temperature of the hole. Let 
f(E u , T) represent the total neutrino flux of energy E v given in eq. 1201 then 
the direct diffuse neutrino flux per unit area today is: 

dN v _ .009 g- 1 [ z ™ [ kT ** d{m) 1 dn 

dE u0 4vr 4 512A J Zmin JkT„(i+z) (/cT*) 5 (1 + zf dz ^ >Z 1 ' 

where E v = E u q(1 + z), z is the redshift at the time of emission, z max corre- 
sponds to max redshift from which particles of energies 100 EeV can reach 
us and z m i n corresponds to minimum redshift. Similarly we calculate the in- 
direct neutrino flux and proton flux from braneworld PBHs by incorpoarting 
fragmentation function in eq. US 



7 Results and Discussion 

The observational constraints on the mass fraction of black holes at evapora- 
tion is given by the quantity a evap = Ppbh( ^ M ^ . All the constraints with reason 

^ Prad 

are given in table 1 of ref. [27|. We have taken upper limits on the PBH 
densities corresponding to their evaporation time for our calculation. The 
results for the neutrino and proton flux from 4d primordial black holes are 
given in Fig. and Fig. |2] respectively. The neutrino flux from the brane 
world PBH is given in Fig. El As we see from Figs. Q and El neutrino flux 
at energy 10 20 eV from PBHs evaporating today is many orders of magni- 
tude smaller compared to the existing neutrino flux limit. For early decaying 
PBHs it can be noticed that neutrino flux at 10 20 eV is even smaller. 

For PBHs of mass M ~ 5 x 10 14 g, neutrino flux at 10 20 eV energy 
is eight orders of magnitude smaller than the neutrino flux limit at that 
energy. If these PBHs are clustered in galactic halos then their density may 
be somewhat higher. For PBHs of mass M ~ 10 13 g and M ~ 5 x 10 10 g, 
neutrino flux at 10 20 eV is very small compared to neutrino flux limit because 
the number of UHE particles emitted by low mass black holes decrease rapidly 
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Figure 3: Direct Neutrino flux today from 5d BWPBHs evaporating at red- 
shift z = O(solid), z = lOOO(dotted) and z = 10 6 (short dashed). Sim- 
ilarly indirect neutrino flux today from 4d PBHs evaporating at redshift 
z = O(dottted), z = 1000 (long dashed) and z = 10 6 (large spaced dots). 

as their mass decreases. Hence, although the constraints on the lower mass 
PBHs are much weaker, the final flux produced by the low mass PBHs turns 
out to be much smaller. For M « 10 13 g the dominant constraint comes from 
entropy production at nucleosynthesis. While for M ~ 5 x 10 10 g the main 
constraint comes from deuterium destruction. 

The astrophysical and cosmological constraints on brane world PBHs are 
obtained in ref. [22 123 121] ■ I n r ef- |2*H1 |2H] constraints on brane world 
PBHs are obtained from high energy diffuse gamma ray and from cosmic ray 
antiproton flux. We observe that neutrino fluxes are much smaller compared 
to 4D PBHs even considering maximum allowed PBH densities at their cor- 
responding evaporation era. This is understandable because the temperature 
of 5D PBHs is small compared to the 4D PBHs of same mass. UHE neutrino 
flux from brane world PBHs is much smaller compared to 4D PBHs of same 
mass. 

We also calculate proton flux from 4D PBHs and brane world 5D PBHs. 
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It turns out that proton flux at 10 20 eV from 4D PBHs decaying today is 
roughly ten orders of magnitude smaller than the UHECR flux. Proton flux 
from early decaying 4D PBHs and brane world PBHs are even smaller. 

The main uncertainities in our calculation include: (a) lack of information 
about the degrees of freedom available at high energies when black hole 
temperature exceeds the energies currently attained in collider experiments 
and (b) extrapolation of the fragmentation function to high energies. These 
unknowns at high energy might change the above picture somewhat but is 
unlikely to change our results qualitatively. 

To conclude, we find that PBHs, decaying in the early universe, contribute 
negligibly to the ultra high energy cosmic ray flux. However it may be 
interesting to repeat our calculations for other superheavy particles, decaying 
in the early universe, which may give a larger contribution. 





10 17 


10 12 


10 4 




6.214 x 1(T 14 


1.338 x 10" 15 


2.885 x 10" 18 


^maxi^CITl) 


.01 


.01 


.01 


M min (g) 


1.549 x 10 9 


4.898 x 10 6 # 


4.898 x 10 2 


M max (g) 


3.175 x 10 14 


1.339 x 10 13 


2.885 x 10 10 


T ■ 

J- mm 


2.3 x 10~ 23 T 4 


2.45 x 10~ 22 T 4 


2.45 x 10- 20 T 4 


T 

J- max 


1.71 x 10~ 2U T 4 


4.053 x 10" iy T 4 


1.88 x 10~ i6 T 4 



Table 1: Mass, temperature ranges for BW PBHs at three different epochs. 



8 Appendix 

Here we give expressions for f(E u ,T) for different cases. 

1. The indirect neutrino flux due to decay of muons and pions may be 
expressed as, 

m,T) = B J f(E,,T) dn ^£ K) dE, , (32) 

where B is the number of degrees of freedom. For the case of muon decay, 
we can express, 

dn v {E^E v ) _ 2 



dE v jm^ 



f(x) (33) 
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where x = 2^ 

f(x) = 2x 2 (3-2x) (34) 

for and 

f(x) = 12x 2 (l -x) (35) 

for z/ e . In the present case B = 4 because we have contributions from fi + and 
Similar calculation can be done for spectrum from pion decay. 
2. For the fragmentation of quarks— *> pions ^muons — ^neutrinos, we have 

/(£„, T) = B, J f(E q , T) d ^§ E -> """^ ^>/»^- *± IEW4E, 

dn " (36) 

where is the neutrino spectrum in a decay and is the muon spectrum 
in a 7r decay. B q is degrees of freedom of all six quarks and six antiquarks. 
Each quark has two degrees of freedom. Hence B q =24 in this case. We can 
express, 

dn u {E m u,E u ) dn u (E^,E u ) dy dy 

dK = dy dE» = {9o{y) ~ P9l{y)) dW» (37) 

in the limit ^> , m e , where 

go(y) = l~3y 2 + \y* (38) 

9i{y) = \-^>y 2 + \y 3 (39) 

for and 

g (y) = 2-6y 2 + Ay 3 (40) 

gi (y) = -2-12y-18y 2 + 8y 3 (41) 

for v e . In eq. E7|P is the projection of the muon spin in the muon rest frame 
along the direction of the muon velocity in the laboratory frame 

P= f -i±I (42) 

where r = • Similarly muon spectrum from pion decay is 

dn^EnEj 1 1 



d£„ 1-r 



(43) 



3. Similar expressions can be obtained in case of neutrinos from W-boson 
decays. 
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